SLAC-PUB-12930 
October 30, 2007 



Relativistic dark matter at the galactic center 



Mustafa A. AminQ 

Dept. of Physics and Kavli Institute for Particle Astrophysics and Cosmology 
Stanford University, Stanford, CA 94305-4060 

and 

TOMMER WlZANSK\FF 

Stanford Linear Accelerator Center 
2575 Sand Hill Road, Menlo Park, California, 94025, USA 

Abstract 

In a large region of the supersymmetry parameter space, the annihilation cross 
section for neutralino dark matter is strongly dependent on the relative velocity of the 
incoming particles. We explore the consequences of this velocity dependence in the 
context of indirect detection of dark matter from the galactic center. We find that the 
increase in the annihilation cross section at high velocities leads to a flattening of the 
halo density profile near the galactic center and an enhancement of the annihilation 
signal. 
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1 Introduction 



Indirect detection is one of the most promising avenues for the discovery of dark matter 
through its non-gravitational effects. Many efforts are underway and more are planned to 
detect the products of dark matter annihilations [HE1EII11E1E1E1IHIE1IID]. The best 
places to look for the annihilation signal are regions where the density of dark matter is 
expected to be high, for example, centers of dark matter halos [H] . center of stars [T2J [13] 
and neighborhoods of compact objects |14j . 

We concentrate on the sub-parsec region around the super-massive black hole (SBH) at 
center of our galaxy (M bh k4x 10 6 M q [TS], [TB]). Gondolo and Silk [T7] argued that a 
sharp dark matter spike should form around the SBH leading to a large enhancement of the 
annihilation signal. Subsequent authors (for example jT8j [19l ED]) qualified this statement, 
pointing out several phenomena which would have the effect of smoothing and reducing the 
spike. The debate over the existence of a dark matter spike at the center of the galaxy has 
yet to be resolved. For the purpose of this paper, we assume that a spike does exist. 

In this study we discuss a new correction to the predictions for the annihilation rate and 
halo profile around the SBH. We point out that near the black hole the dark matter particles 
will be moving sub-relativistically (v/c < 0.2). This is in contrast to the usual assumption 
whereby the dark matter is taken to be cold and slow. In fact, most previous calculations 
(see for example [T7]) have been performed in the limit (v/c) — > where v is the relative 
velocity between particles. For a certain class of supersymmetric dark matter models, the 
cross section for annihilation can be enhanced by several orders of magnitude in the vicinity 
of the SBH due a strong dependence on v. In the presence of a central dark matter spike this 
can produce a measurable correction to the observed annihilation signal. In addition, the 
enhanced cross section leads to depletion of the spike and a widening of the "annihilation 
core". We explore these two effects for a variety of spike profiles to account for the many 
astrophysical uncertainties regarding the nature of the density profile. 

We find that the enhancements in the annihilation signals occur primarily in models 
for which the indirect detection signals are too small to be seen by current experiments. 
However, these models are quite plausible theoretically and are even preferred by some 
criteria. We can easily imagine a scenario in which particle physics experiments point to one 
of these theories as a correct description of nature. This will motivate dedicated gamma ray 
observations concentrating on objects where dark matter is likely to be concentrated. We 
will argue that, in this situation, the velocity-dependent enhancement of the annihilation 
cross section must be taken into account. 

The rest of the paper is organized as follows. In Section [2] we give a brief review of 
supersymmetric dark matter and enumerate the circumstances whereby a strong enhance- 
ment to the annihilation cross section may arise. In Section [3] we estimate the corrections 
to the halo profile arising from the enhanced annihilation rate and calculate corrections to 
the annihilation signal. Our conclusions are presented in Section |4| 
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Figure 1: A typical set of Feynman diagrams contributing to the self annihilations of a neu- 
tralino into Standard Model particles. The LSP is denoted by xii Xj is a heavier neutralino 
and A is the neutral CP-odd Higgs boson. 

2 Supersymmetric dark matter 

For the purpose of this study we restrict ourselves to the minimal supersymmetric standard 
model (MSSM). In this class of theories there exist four neutral fermionic mass eigenstates - 
the neutralinos. The lightest of these is often the lightest superpartner in the theory (LSP) 
and provides a good dark matter candidate. We are interested in describing the conditions 
under which the annihilation of the LSP to standard model particles exhibits a strong velocity 
dependence leading to an enhancement of the indirect detection signal. 

A sample of the most important Feynman diagrams contributing to neutralino anni- 
hilations are depicted in Figure [TJ First, a pair of neutralinos may exchange a fermion 
superpartner (sfermion), producing two standard model fermions. Fermions may also be 
produced through an s-channel exchange of a heavy scalar, in this case the A Higgs boson. 
Notably, this diagram does not admit a p-wave component, a fact which will be important 
in the coming analysis. Finally, the neutralinos may annihilate to standard model gauge 
bosons. In Figure [I] we present the annihilation to two Z° bosons via the exchange of a 
heavier neutralino. 

In the MSSM, neutralinos are Majorana particles. This leads to a well-known helicity 
suppression of the amplitude for pair annihilation into light fermions [21] - If X denotes the 
dark matter particle and / the fermion, the s-wave cross section will be suppressed by a 
factor of 

2 2/2 

z = mj/m x . 

For the annihilation of a 200 GeV neutralino to Standard Model leptons, z 2 is less than 10~ 4 . 
Consequently, the p-wave annihiliation, which is suppressed only by v 2 /c 2 , may dominate. 
We can therefore conclude that for models where the LSP annihilates primarily to fermions, 
the annihilation cross section will exhibit a strong velocity dependence. It is this effect 
which lies at the core of our present work. In the next few paragraphs, we review how these 
helicity-suppressed cross sections arise in the MSSM parameter space. A more complete 
description can be found, for example, in [22]. For an excellent review of the MSSM and 
supersymmetry in general see [23J. 
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Figure 2: <7t> as a function of (v/c) for sample models from the f coannihilation region. In 
(a) rriA is scanned, increasing from top to bottom, and in (b) mf is scanned, increasing from 
top to bottom. 



As mentioned, only annihilations to fermions undergo helicity suppression. We would 
like to identify the regions of parameter space for which the dominant annihilation channels 
do undergo helicity suppression and the resulting process is p-wave. It is this class of models 
which will exhibit a strong velocity dependence. 

In the MSSM, each neutralino is a linear combination of the superpartners of two neutral 
gauge bosons and two neutral Higgs bosons. It is typically parametrized by 

Xl = Z a B° + Z l2 W° + Z & H* + Z l4 H° 2 , 

where Xi is the i th neutralino and tildes denote superpartners. The partners of the B°, W° 
and Higgs bosons are usually called wino, bino and Higgsino respectively. Of these four only 
the bino is a gauge singlet, meaning that it does not interact with gauge bosons. Thus, by 
making Z\\ large compared to the other components, we can eliminate the third diagram in 
Figure [TJ leaving only fermionic processes. The annihilation of Majorana particles through 
a scalar coupling can only take place in the s-wave. Thus, if the second diagram were to 
dominate over the first, the cross section would indeed be helicity suppressed but no strong 
velocity dependence would arise. To suppress this diagram we demand that the A boson 
is significantly heavier than the fermion superpartners. We must also make sure that no 
resonance enhances the A diagram, that is, cannot be too close to 2m x . 

We may now ask how generic are these constraints? The condition of large Z\\ is quite 
generic. The theoretically compelling assumption of gauge unification naturally leads to a 
bino that is lighter than the wino by a factor of two [21]. In many classes of models, for 
example, minimal supergravity, the condition of electroweak symmetry breaking requires 
that the higgsinos are quite heavy. In these cases the LSP is mostly bino. The mass of the 
A boson is, in principal, unconstrained and can easily be large enough to suppress the A 
exchange diagram. 

Thus, helicity-suppressed dark matter annihilation is quite likely in models of supersym- 
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metry. This implies a strong dependence of the annihilation rate on the relative velocity of 
the incoming particles. In the following, we will consider sample MSSM models from this 
region of parameter space. We will choose models for which the predicted dark matter relic 
density agrees with value Q x h 2 w 0.1 given by the WMAP experiment [25]. All cross sections 
and relic density calculations were performed using the DarkSUSY software package [26] . 

In Figure [2j we show the annihilation cross section times the velocity for the dark matter 
particle in several of these models. The strong velocity dependence is evident, enhancing the 
total cross section by several orders of magnitude over the value at v = 0. We show how the 
low and high velocity behaviors of the cross section can be tuned independently by varying 
different supersymmetric parameters. In Figure [2] (a) the mass of the the A is scanned. 
As it is decreased the s-channel diagram in Figure [T] becomes increasingly important, the 
s-wave component of the amplitude increases, and the v = cross section grows. In[2](b) the 
mass of the stau is scanned. As this parameter is decreased the t-channel scalar exchange 
diagram is enhanced and the p-wave, velocity dependent, component of the cross section 
grows. Notably, we must vary the mass of the LSP along with the stau mass in order to 
maintain the correct relic abundance. 

The curves in Figure [2] can be well fit by expressions of the form 



where (<jv)q and (crv)i are fit parameters. We use this approximation in the following sec- 
tions. 



3 Astrophysics 

In the previous section we introduced a class of SUSY models for which the neutralino 
annihilation cross section shows a strong velocity dependence. In this section we explore the 
consequences for the density profile and the expected annihilation signal from the dark matter 
in the sub-parsec region around our Galactic center. To the best of our knowledge, this effect 
has been ignored in the literature. We find that neglecting this velocity dependence leads to 
an underestimation of the size of the annihilation core as well as the expected annihilation 
signal. In addition, we find that these corrections depend strongly on the density profile of 
the halo. 

3.1 Density profile 

The annihilation signal depends on the the density profile of dark matter. In particular, it 
is sensitive to the profile in the sub-parsec region around the central black hole where the 
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density is expected to be quite high. Our first goal is to understand how the density profile 
changes in this region when we include the velocity dependence of the cross section. 

The density profile depends on a number of physical processes such as the initial phase 
space distribution of the dark matter particles that collapsed to form the halo, the steepening 
of the profile due to the baryons, scattering by stars, loss to the central black hole, black 
hole or galactic merger history etc. A detailed calculation of the density profile is beyond 
the scope of this paper (see [27] for an excellent review). We assume the following density 
profile [30] 

(p(r c ) 10r g <r <r c 

Po (r/r hh y Jsp r c < r < r bh , (2) 

Po (r/r hh y a r bh < r , 

where the core radius, r c , is defined in the next paragraph. For r < 10r 9 the density of 
dark matter particles decreases rapidly and vanishes at r = 4r g [IT] . In the above expression 
r g w 4 x 10~ 7 pc is the Scwarzchild radius of the central black hole and r b h ~ 2pc is the 
radius at which the mass of the stars within that radius is twice the mass of the central black 
hole. We take p(r b h) = Po ~ 100M o pc~ 3 , though it could be higher [30] ■ We parametrize 
our ignorance regarding the nature of the profile using the two coefficients a and 7 sp . 

We now turn to the core radius r c . As discussed above, the density profile is determined 
by self annihilation, scattering by stars, loss to the SBH etc. Scattering by stars drives the 
density profile to a power law [27]. If the density gets too high, annihilation becomes efficient 
enough to prevent further rise in the density. This results in the formation of a flattened 
core near the galactic center. The radius at which the core starts forming is determined by 

Tannic) ~ t heat (3) 

where t bea t ~ 2 x 10 9 yrs [29j [27J is the timescale for heating of the dark matter particles due to 
scattering by stars. The annihilation rate r ann (r) = p(r)av(r)/m where m is the mass of the 
dark matter particle. The position dependence of av arises due to its velocity dependence. 
For a virialized halo, (v/c) 2 ~ r g /2r). Since the dark matter density is significant for 
r > I0r g , the relavent velocities are bounded by iy/c) 2 < 0.05. 

We consider a model taken from the stau coannihilation region of mSugra. The mass of 
the LSP m = 166 GeV and the mass of the lightest stau is 173 GeV. The relic density is 
O dm /j 2 « 0.1. In this model, (av) = 9 x 10 _30 cm 3 s _1 and (ca>)i = 8.9 x 10~ 26 cm 3 s _1 (see 
equation ([I])). We will refer to this model as our fiducial model. Whenever a parameter is 
not explicitly defined or varied, its value is taken from this model. 

For our fiducial model with 7 sp = 3/2 and a = 1, the core radius is r c m 14r s , with 
a core density of p(r c ) ps 2 x 10 8 p - If we ignore the velocity dependence, then we do not 
get a core. In Figure |3](a) we plot r c for different (av) and (av ) 1 for the same density 
profile. We note that the size of the core is not independent of (cv)\ and ignoring it leads 
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Figure 3: (a) The variation of the core radius with (av)o and (av)\. Ignoring the (av)\ leads 
to an underestimate of the core radius, (b) Change in the spike profile as function of 7 sp . 
Note that for large values of 7 sp , the radius of the annihilation core is also large. The dotted 
lines indicate the density profile when we set (av )i = 0. 

to an underestimation of the core size. If {<jv)i has a dominant contribution in determining 
the core radius r c , it has to be significantly larger than (<jv)q. This is because the factor 
(v/c) 2 ~ r g /2r in front of (crv)i is small unless we are close to the central black hole. 

Another important factor that determines the size and density of the core is the steepness 
of the density profile parametrized by 7 sp . For the fiducial model, the size of the core as 
well as the density increases with increasing 7 sp as shown in Figure |3](b). The dotted lines 
represent the density profiles for (av )i = 0. 



3.2 Annihilation Flux 



The flux of photons (as observed by us) can be written as 

> 4 f^': 1 <4) 

2m 2 J 4vr|d + r| 2 v ; 

where p(r) is the dark matter density and av(r) is the annihilation cross section times the 
typical relative velocity of the annihilating particles. In the above expression d is the vector 
joining the sun and the galactic center, m is the mass of the annihilating dark matter particles 
and N is the number of photons (above the detector thresh-hold) produced in the annihilation 
process. The integral is done over a solid angle which depends on the angular resolution of 
the detector. We take this to be AQ « 10~ 5 sr which is the approximate angular resolution 
for GLAST. We remind the reader that the position dependence of Nav (r) arises from the 
position dependence of the velocity in a virialized halo. Due to this position dependence of 
the cross section, we cannot simply separate the particle physics and astrophysics aspects of 
the integral as is commonly done in the literature. 
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Figure 4: (a) The annihilation flux as a function of (crv)o and (crv)i. $ is expressed in units 
of 6 x 10 _10 cm _2 s _1 which is the (approximate) sensitivity of GLAST at an energy threshold 
of 1 GeV. (b) The ratio of annihilation fluxes $/$o where '0' refers to the flux calculated by 
setting (crv)x = 0. 



The annihilation signal depends on the cross section in two ways: Explicitly through 
Nav appearing in equation Q and implicity through p which depends on av as discussed 
in the Section 13.11 

For the fiducial model with 7 sp = 3/2 and a — 1, the expected flux is enhanced by a 
factor of ~ 4 compared to the case when the velocity dependence is ignored. One should 
view this number with caution, since it does depend strongly on the parameters of the model. 

In Figure |4]^a), we plot the annihilation flux as a function of (<jv)q and (av)\ with the same 
halo profile. In Figure |4^b) we plot the ratio of the fluxes, with and without the velocity 
dependence in the cross section : $/$o where '0' indicates that we set (av)\ = 0. As 
expected, ignoring the velocity dependence of the cross section leads to an underestimation 
of the flux. The enhancement is large when (av)i/(<Jv)o is large. 



Next, in Figure [5] we show $/$o as a function of 7 sp . For 7 sp < 1.6 the enhancement 
increases with 7 sp . This is due to the increase in the density very close to the blackhole 
compared to regions farther away which results in a greater fraction of dark matter particles 



having high velocities. However, as we saw in Section 3.1 the size of the annihilation core 
also increases with 7 sp . When the core becomes sufficiently large (7 sp > 1.6), the spike profile 
is significantly flattened (see Figure 3(a)). This decreases the enhancement of the flux. 

We note that the enhancement of the signal occurs in models that are not detectable by 
current or planned experiments. For our fiducial model, the flux is two orders of magnitude 
below GLAST sensitivity (see for example [31]). This is mainly due to the small (crv)o 
since it is (crt>) that determines the annihilation flux in regions with r > 10 4 r 9 . It is 
tempting to explore the SUSY parameter space with the aim of finding models with a large 
(crv)o and (av )i/(crv )o, so that the flux is large to begin with and the velocity dependent 
enhancement provides a further boost. However, relic dark matter abundance constrains 
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Figure 5: (a) Flux as a function of 7 sp (b) The ratio of the annihilation fluxes: $/$o a s 
a function of 7 sp . '0' indicates that we set (av)i = 0. The enhancement decreases after 
7s P ~ 1.6 due to a significant increase in the size of the annihilation core. 

(av)i(v /c) 2 < 10~ 26 cni 3 s~ 1 . Thus, for (av)i/(av) > 10 4 , (crv) is typically small leading to 
a small overall flux. 



4 Discussion 

In this paper we have discussed the consequences of relativistic dark matter near the black 
hole at the center of our galaxy. We have argued that, in general, the commonly used 
approximation whereby the relative velocity of dark matter particles is taken to vanish may 
be inapropriate. In regions very close to the black hole, the cold dark matter is no longer 
cold. If the dark matter has accumulated in a sharp spike around the black hole, this 
region may account for a large fraction of the expected signal. We presented a specific class 
of supersymmetric models in which the dark matter annihilation cross section is strongly 
dependent on the relative velocity of the incoming particles. In these theories, the expression 
for the annihilation flux no longer separates neatly into factors depending on the astrophysics 
and the particle physics. When the full velocity dependent cross section is considered, the 
annihilation flux receives up to an order of magnitude enhancement over the v = value. 
In addition, we found that the enhanced cross section effects the halo profile close to the 
galactic center. The increased annihilations deplete the spike and widen the annihilation 
core. 

We explored the the change in the density profile and annihilation signal for annihilation 
cross sections of the form av = (crv)o + (<jv)i(v/c) 2 . We showed how the annihilation core 
size and the flux changed as a function of (av) and (crv)i. To account for the astrophysical 
uncertainties in determining the dark matter density near the galactic centre, we presented 
our results for a variety of spike profiles. 
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None of the models we have considered are detectable by current or upcoming gamma ray 
observations. If the neutralino is the dominant component of dark matter and is produced 
thermally, the cross section at high velocity cannot be larger than about 10 -26 cm 3 /sec; 
otherwise the relic abundance would be too small. In most regions of the galaxy today, the 
neutralino velocity v/c would be very small, and the annihilation signal would be highly 
suppressed. However, if particle physics observations should indicate a scenario like those 
we have described, it would be worthwhile to mount dedicated gamma ray observations 
concentrating on the galactic center and the centers of nearby galaxies. Uniquely in those 
environments, in the neighborhood of the central black holes, the annihilation cross section 
would be enhance by the effect described in this paper. 
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6 Appendix 

In this appendix we provide some analytic approximations to the flux integral, equation Q. 
We will assume that the velocity dependence of the cross section takes the form of equation 
(Jl]), although this is not essential in the numerical calculations. 

We split the flux integral into three parts; $ = $ CO re + $s P ike + $haio bases on the density 
profile ([2]). In most cases, the largest contribution to the signal comes from the spike. 
However, the contribution from the core and halo is not always negligible. For example in 
the fiducial model, the spike, core and halo contribute 74, 15 and 1 percent of the signal 
respectively for an angular resolution of AQ = 10~ 5 sr and 7 sp = 3/2, a = 1. 

For the density profile, equation ([2]), we can calculate the core and spike parts of the 
integral analytically (since r c ,r bh <C d). For 7 sp = 3/2, the flux from the core and spike is 
given by (the halo integral is harder to do analytically, unless one assumes that most of the 
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flux comes from r -C d) 

_ 1 (Nav)o p 2 Q rl h 



spike 



6 d? 

j (Nav)o p 2 rl h ^ ( r bh 



2d 2 m 2 




(5) 



It is important to note that r c depends on the cross section. For the case when 7 sp ^ 3/2 
we have 



l(iVai;)oP^ c 3 



6 d 2 m 2 



10rv 



3(iW)i?~ g 



10r c 



4 (M™) r c 
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rbh) 
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3 - 27 sp (Nav) 1 r g 



''bh 
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4(1 - 7sp ) (Nav)o r bh 



Again r c and p c depend on the cross section and 7 sp . 
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